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Abstract 

We present the calculation of the atmospheric neutrino fluxes for the neutrino experiments 
proposed at INO, South Pole and Pyhasalmi. Neutrino fluxes have been obtained using ATMNC, 
a simulation code for cosmic ray in the atmosphere. Even using the same primary flux model and 
the interaction model, the calculated atmospheric neutrino fluxes are different for the different sites 
due to the geomagnetic field. The prediction of these fluxes in the present paper would be quite 
useful in the experimental analysis. 
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I. INTRODUCTION 



The study of neutrino physics with atmospheric neutrinos has a long history with the 
first observations of muons produced by in deep underground laboratories of Kolar Gold 
Field(KGF) mines in India jl| and East Rand Propietary Mines (ERPM) in South Africa J2I. 
It was the Kamiokande 3, 1MB fl and some other atmospheric neutrino experiments [5| 
which gave a clear evidence of a deficit in the atmospheric muon neutrino flux. With the 
Super-Kamiokande [6( experiment, it is now well established that atmospheric neutrinos do 
oscillate. Both theoretically and experimentally lots of activities are going on to precisely 
determine neutrino oscillation parameters and experiments are planned with atmospheric as 
well as with the accelerator, reactor and solar neutrinos. 

India is going to establish a neutrino observatory at the Theni district of Tamilnadu Q. 
This India-based Neutrino Observatory (INO) is proposed to be installed at 9°96'7"N, 
77°26 7 E, and their primary aim is to study atmospheric neutrinos in a 1300 meters deep 
cave. The detector that is planned to be used is a magnetized iron calorimeters (ICAL) 
of 50kT. At the South Pole, PINGU (Precision IceCube Next Generation Upgrade) is pro- 
posed as the upgrade version of IceCube (-90°, 0°) jg] to observe atmospheric neutrino at 
lower energies by increasing the optical modules many times of the present setup, under 
the surface of the ice. Another new generation underground neutrino experiment facility 
called LAGUNA(Large Apparatus studying Grand Unification and Neutrino Astrophysics) 
is proposed at Pyhasalmi Mine [9| in Finland, which is among the deepest mines in Europe 
at 63°39'N, 26°02'E. 

In this paper, we have calculated atmospheric neutrino fluxes for these three sites in a 
3D scheme using ATMNC (ATmospheric Muon Neutrino Calculation) code 10h12| for the 
cosmic ray propagation in atmosphere with JAM, which is used in PHITS (Particle and 
Heavy-Ion Transport code System) 13j, in the hadronic interaction at lower energies (< 
32 GeV) . The JAM interaction model agrees with the HARP experiment [l^] a little better 
than DPMJET-III [15]. Earlier ATMNC code has been applied to the study of muon flux at 
several altitudes, at sea level, mountain altitude, and at balloon altitudes, where accurate 
measurements exist. The Monte Carlo generator with JAM shows a better agreement than 
the former one without JAM. Then it is applied to the calculation of atmospheric neutrino 
flux for several sites as the GranSasso, SNO, Kamioka and others 10h12[. 
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In spite of using the same primary flux model and the interaction model for the different 
sites, the calculated atmospheric neutrino fluxes are different due to the geomagnetic field. 
The geomagnetic field affects cosmic rays both inside and outside of the atmosphere. First, 
it acts as a filter for low energy cosmic rays, and secondly, it deflects the charged particles 
in the atmosphere. These two effects are mainly controlled by the horizontal component of 
the geomagnetic field. It would be interesting to study the atmospheric neutrino flux for the 
three sites with different position in the geomagnetic field. In Fig. [H we show the strength 
of the horizontal component of geomagnetic field obtained using IGRF2010 model 16| with 
the position of these sites. It can be seen that the INO site is close to the region where 
the strength of the horizontal component of geomagnetic field is the largest on the earth. 
The South Pole is close to the magnetic pole, where the horizontal component is zero. At 
Pyhasalmi mine, the horizontal component of geomagnetic field is also small but not zero 
as it is little far away from the magnetic pole. 

We proceed as follows. In section-HTl we present the main features of the calculation 
scheme and in section 4ITTl the results of the atmospheric neutrino fluxes have been shown 
and discussed. Finally, in section-HVt we summarize the results and conclude our findings. 



II. CALCULATION SCHEME 



The scheme for calculating the atmospheric neutrino fluxes has been discussed in detail 
in the earlier work 
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model based on AMS 
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121 ] . We present here the main features. We use the primary flux 



18| and BESS 
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20J data. We use IGRF2010 (l6| model for 



the geomagnetic field. For the atmosphere model, we use NRLMSISE-00 21[ instead of US- 
standard76 [22j which was used in our earlier works 10Hl2|. Actual calculation is carried 
out in the Cartesian coordinate system which has the origin at the center of the Earth, 
with the Z-axis extending to the north pole, and we consider the surface of the Earth to be 
a sphere with a radius of R e = 6378.180 km. However, the position on the Earth is well 
represented by the spherical polar coordinate system (r, 9, <f) with r = R e which is related 
to the Cartesian coordinate system by 



x 



= R e sin 9 cos0, x = R e sin 6* sin</> and z = R e cos 9. 



The local coordinate system at the detector is constructed based on this polar coordinate 
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system. The direction of the x-axis is in the increasing direction of 9, the direction of the 
y-axis is in the increasing direction of <f), and the direction of the z-axis is in the increasing 
direction of r. Therefore, the azimuth angle is measured counterclockwise from south in 
the local coordinate system. In addition to the surface of the Earth, we assume three more 
spheres; the injection sphere, the simulation sphere, and the escape sphere. We have taken 
the radius of the injection sphere as R\ n j = R e + 100km, and the radius of simulation 
sphere and the escape sphere are taken to be R esc = R S i m = 10 x f? e =63781.80 km. 

Cosmic rays are sampled on the injection sphere uniformly towards the inward direction, 
following the given primary cosmic ray spectra. Before they are fed to the simulation code 
for the propagation in air, they are tested to determine whether they pass the rigidity cutoff 
or not. For a sampled cosmic ray, the "history" is examined by solving the equation of 
motion in the negative time direction. When the cosmic ray reaches the escape sphere 
without touching the injection sphere again in the inverse direction of time, the cosmic ray 
can pass through the magnetic barrier following its trajectory in the normal direction of 
time. The propagation of cosmic rays is simulated in the space between the surface of Earth 
and the simulation sphere. 

We use the JAM interaction model for hadronic interactions below 32GeV, as this shows 
a better „t with the HARP expert Q and it agrees with the observed m uo„ 
flux at sea level, at mountain altitudes and at balloon altitudes. For energies above 32GeV, 
we use DPMJET-III 15[ interaction model. We have checked the smooth interpolation when 
switching from the JAM model to the DPMJET-III interaction model. 



III. RESULTS AND DISCUSSIONS 



In this section, we present the results of simulation for the atmospheric neutrinos at INO, 
South Pole and Pyhasalmi sites. First we present the results for the atmospheric neutrino 
fluxes as a function of azimuthal angle (p. These results are presented for u^, z> M , v e and v e 
for (anti)neutrinos of 3.2GeV. In Fig. [2j we present the results for the INO site, in Figj3] 
the results are presented for the South Pole site and in Fig. H] the results are presented for 
the Pyhasalmi site at -E J/ =3.2GeV. In these figures we show the variation of atmospheric 
neutrino flux as the function of the azimuthal angle averaging them over the five zenith angle 
ranges,! > cos6> > 0.6,0.6 > cos# > 0.2, 0.2 > cos# > -0.2,-0.2 > cos# > - 
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FIG. 1: Magnitude of the horizontal component of geomagnetic field in IGRF2010 model fit] ]. 
Square stands for the position of India-based Neutrino Observation (INO) site, diamond for the 
Pyhasalmi mines, and bottom bar for the South Pole. 

0.6 and -0.6 > cos 9 > -1. We find that the variation of the atmospheric neutrino 
flux has a complex structure at low (anti) neutrino energies, due to the rigidity cutoff and 
muon bending in the geomagnetic field. This variation remains almost the same for the near 
horizontal direction even above 10 GeV. Due to the high rigidity cutoff at the INO site this 
variation is more complex than the other two sites discussed here and for the South Pole 
site this variation is the least. 

Furthermore, the variation of upward going neutrinos is much more complicated than 
the variation of downward going neutrinos. This is due to the fact that the upward-going 
neutrinos are produced in wide area on the Earth, and there are large variation of rigidity 
cutoff and geomagnetic field. On the other hand, the downward going neutrinos are produced 
just above the detector. 

In Figs. El |6]and [71 we present the results for the atmospheric neutrino fluxes as a function 
of the zenith angle after averaging over all the azimuthal angles. The results are presented 
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FIG. 2: The azimuthal angle dependence of atmospheric neutrino flux, averaged over zenith angle 
bins of 1 > cosr? > 0.6 (solid line), 0.6 > cos 9 > 0.2 (long dashed), 0.2 > cos 6 >-0.2 
(short dashed), -0.2 > cos# > -0.6 (dashed-dotted), and -0.6 > cos 9 > -1 (dashed double- 
dotted), calculated for the INO site at (anti) neutrino energy E=3.2GeV. 



for neutrino energy of lGeV, 3.2GeV and lOGeV respectively, and for the three sites viz. 
INO [7], South Pole [8j and Pyhasalmi {9] in each of these figures. At 1 GeV, there are large 
up-down asymmetries in the atmospheric neutrino flux at all the three sites. The downward 
going neutrino flux is larger at the South pole Q| and Pyhasalmi |9f sites, while upward 
going neutrino flux is larger at the INO site [7] due to the different rigidity cutoff. These 
asymmetries decrease with the increase in neutrino energy, and almost disappear at 10 GeV. 
However, there appears an up-down asymmetry at South Pole, due to the difference in the 
observation altitude. 

We note that the horizontal/vertical flux ratio for the INO site is much smaller than 
for the other sites even at 3.2 GeV. This could be understood by the fact that the rigidity 
cutoff still affect the neutrino flux at 3.2 GeV at the INO site, and the rigidity cutoff is more 
effective in the horizontal direction. 
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FIG. 3: Same as Fig. [2] at (anti)neutrino energy E=3.2GeV for the South Pole site. 
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FIG. 4: Same as Fig. [2] at (anti)neutrino energy E=3.2GeV for the Pyhasalmi site. 
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FIG. 5: The zenith angle dependence of atmospheric neutrino flux at E=lGeV, averaged over 
all azimuthal angles calculated for INO, South Pole and Pyhasalmi sites. Here 9 is the arrival 
direction of the neutrino, with cos 9 = 1 for vertically downward going neutrinos, and cos 9 = — 1 
for vertically upward going neutrinos. 

In Fig. [8], we have shown proton spectra which produce neutrino above 3.2 GeV, from the 
near vertical (1 > cos9 zenith > 0.4) and near horizontal (0.2 > cos 9 zenit h > —0.2) directions 
for the INO site, with effectively no rigidity cutoff site (SNO [23[) and intermediate rigidity 
cutoff site (SK [6]). On comparing the proton spectra for these sites, it may be noticed that 
the rigidity cutoff works more efficiently for the near horizontal direction than for the near 



vertical direction, especially for the INO site [7 



neutrino at South pole 
site y. 



and Pyhasalimi sites 



. The rigidity cutoff for downward going 
9| are effectively the same as that of SNO 



In Fig. |9j we present the results for the atmospheric neutrino spectra averaged over zenith 



a, 



and azimuth angles, for (anti)neutrino energies from 0.5 GeV to fO GeV, for the INO 
South Pole [§] and Pyhasalmi [9] sites. We find that when the neutrino flux is integrated 
over all the angles, the difference in the flux at the South Pole and Pyhasalmi sites which 
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FIG. 6: Same as FigSat E=3.2GeV. 




/ . ■ ■ 



', \ 
'. \ 



0.1 



0.01 



I 1 I 1 I 1 q l L 1 | 1 I 1 I 1 




South Pole 





0.1 




Pyhasalmi 



0.01 1 1 1 L 



I I I 



0.5 1 -1 -0.5 0.5 1 -1 -0.5 0.5 1 
cos9 cos6 cos9 



FIG. 7: Same as FigEJat E=10GeV. 
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FIG. 8: The spectra of protons which create the neutrino above 3.2 GeV from the near vertical (1 > 
cos 6 zenith > 0.4, left panel) and near horizontal (0.2 > cos 9 ze nith > —0.2, right panel) directions 
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23l|). intermediate rigidity cutoff site (SuperK [6]) and 




E (GeV) 



FIG. 9: Atmospheric neutrino flux sum over all zenith and azimuthal angles for the INO, 
South Pole and Pyhasalmi sites. These results are summed over all the flavors of neutrinos viz. 



+ Vu + V e + V e . 



appeared in Figs.[5l|6]and[71 disappears. Due to the strong effect of the horizontal component 
of the geomagnetic field atmospheric neutrino flux at the INO site is almost 30% smaller 
at lGeV, than the flux at the two other sites discussed here, while this difference becomes 
smaller with the increase in the neutrino energy, for example, at 3GeV, this difference is 
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10%, but the flux reduces almost by a factor of 20 as compared to the flux at lGeV. 

IV. SUMMARY AND CONCLUSION 

We have studied the atmospheric neutrino flux for the INO, South Pole and Pyhasalmi 
sites using ATMNC with JAM interaction code below 32 GeV. We find that the atmospheric 
neutrino flux is quite different in nature particularly at low and intermediate energies, de- 
pending on the position in the geomagnetic field, and the strength of the horizontal compo- 
nent of the geomagnetic field is a good measure of the deviation in the fluxes at the different 
sites. Especially, the difference of the zenith angle dependence of the atmospheric neutrino 
flux is important in the analysis of the neutrino oscillation. The difference is large at lower 
energies up to the neutrino energies of a few GeV. These results would be useful in the 
analysis of the atmospheric neutrino experiments proposed at these sites. 
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